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Abstract 

The study of nuclear effects for J/-!/; production in proton- nucleus collisions is crucial for a correct 
interpretation of the J/V' suppression patterns experimentally observed in heavy-ion collisions. 
By means of three representative sets of nuclear parton distributions, the energy loss effect in 
initial state and nuclear absorption effect in final state are taken into account in the uniform 
framework of the Glauber model. A leading order phenomenological analysis is performed on J/^ 
production cross section ratios Riy/Bei^F) for the E866 experimental data. The J/ip suppression 
is investigated quantitatively due to the different nuclear effects. It is shown that the energy loss 
effect with resulting in the suppression on Riy^^^lxp) is more important than the nuclear effects 
on parton distributions in high xp region. The E866 data in the small xp keep out the nuclear 
gluon distribution with a large anti-shadowing effect. However, the new HERA-B measurement 
is not in support of the anti-shadowing effect in nuclear gluon distribution. It is found that the 

abs 

xp in the region xp > 0.2. 



J/tp-imcleon inelastic cross section a'^[f depends on the kinematical variable xp, and increases as 
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I. INTRODUCTION 



The production of J/^l) in high energy coUisions has attracted the extensive attention 
from both the nuclear and particle physics communities'^!. J/ if) suppression is considered 
as a most rehable signature for the formation of Quark Gluon Plasma (QGP) in heavy 
ion collisions during the present time. However, there are many different nuclear effects 
in nucleus-nucleus collisions. A robust interpretation of the experimental data in heavy ion 
collisions requires to understand deeply and quantitatively the basic mechanisms responsible 
for the suppression of J/V' production due to the nuclear effects. Therefore, it is desirable 
that a good basehne should be established by means of the study on J/'^ production in 
proton-nucleus collisions to clarify the conventional nuclear suppression mechanism. 

Several proton-induced fixed target experiments (such as, NAS^l, E772[^l, £866^, NASQI^I 
and HEAR-BI^I) have studied the nuclear dependence of J /iIj production cross-sections as a 
function of the Feynman variable xp- These experimental measurements employ the usual 
parametrization by a power law: UpA — (7pN ■ where is the proton-nucleon cross 
section and UpA is corresponding proton-nucleus cross section for a target of atomic mass 
number A. The E866 collaboration W pubhshed the precise measurement of the suppression 
factor a for 800 GeV protons incident on iron and tungsten nuclear targets, relative to 
beryllium nuclear targets with very broad coverage —0.1 < xp < 0.95. The observed 
suppression is smallest at vahies of 0.25 and below, and increases at larger values of 
xp- Recently, the HERA-B collaboration reported the first measurement of nuclear effects 
in J/%1) production extending into the larger negative part of the Feynman-x spectrum, 
—0.34 < xp < 0.14 by 920GeV protons incident on carbon, titanium and tungsten targets. 
The E866 and HERA-B results are compatible within statistical and systematic uncertainties 
in the overlap region, and are systematically above the NA50 results'^! with based on lower 
energy collisions. At lower values oixp, the HERA-B a{xp) measurement indicates a reversal 
of the suppression trend seen at high xp: the strong suppression established by previous 
measurements at high xp turns into a slight tendency towards enhancement in the negative 
Xp region. 

It is now clear that J/ip production occurs in two different steps where a charm quark 
pair is produced first through the interaction of a projectile on a target parton, followed by 
the non-perturbative formation of the colorless asymptotic state. As a consequence, color 
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octet as well as singlet cc states contribute to production. Two formalisms have been 
proposed to incorporate these features, which are the non- relativistic QCD (NRQCD) 1^1, and 
the color evaporation model (CEM) [^1 , respectively. The CEM shows some similarities with 
NRQCD. In CEM calculations, J /ip production via color-octet processes is allowed, hence the 
kinematic dependence of the cross section is similar to NRQCD in the production of J /il^. 
The CEM and NRQCD predictions have been successful in charmonium phenomenology. 
Therefore, the CEM and NRQCD are used to date as a model for charmonium production. 
It is interesting to notice that the measured J/'4> polarization favors the CEM and is not 
within the predictions of NRQCD, and the CEM has fewer free parameters than NRQCD'^l 

The suppression of J /tij production in proton-nucleus collisions could be caused by many 
different nuclear effects. The nuclear effects are usually subdivided in "initial state effects" 
related to the projectiles which produce a charm quark pair and "final state effects" from 
interactions of the charm quark pair or the fully formation oi a J/ip in the nuclear envi- 
ronment. The nuclear effects on parton distribution functions and initial state energy loss 
effect are considered as the most important initial state effects, while the so-called nuclear 
absorption of the charm quark pair is the main final state effect. 

After the production of a charm quark pair, the prc-mcson charm quark pair will travel 
through the nuclear environment. They are subject to the strong interactions with the 
nuclear matter. Such interactions, which reduce the probability that the pre-meson state will 
ultimately form a, J/t/j without being absorbed, can be described by the Glauber model'^^l 
The nuclear absorption cross section is introduced to express the J/ ip production in proton- 
nucleus collisions. 

After the discovery of the EMC effect^^^^, it is well known that nuclear parton distribution 
functions are different from those in a free nucleon. The nuclear modifications relative to the 
nucleon parton distribution functions, are usually referred to as the nuclear effects on the 
parton distribution functions, which include nuclear shadowing, anti-shadowing, EMC effect 
and Fermi motion effect in different regions of parton momentum fraction. As for the nuclear 
parton distributions, the global analysis method has been proposed in the recent years. So 
far, three groups have presented global analysis of the nuclear parton distribution functions 
analogous to those of the free proton. These are the ones by Eskola et al. ( EKS98 ^'^'^'^ and 
EPS[^3])^ by jjjj^^j (HKMl^^l, HKN04115] and HKNOTl^^l), and by de Florian and Sassot 

(nDS) [^^1 . It is noticeable that the nuclear valence quark distributions from different global 
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analysis are nicely in good agreement. The nuclear sea quark and gluon distributions in 
general are still quite badly constrained, resulting in large differences between the different 
sets. 

The initial state energy loss effect is another nuclear effect apart from the nuclear effects 
on the parton distribution as in deep inelastic scattering. In high energy proton-nucleus 
scattering, the projectile rarely retains a major fraction of its momentum in traversing the 
nucleus. The projectile can lose a finite fraction of its energy due to the multiple collisions 
and repeated energy loss in the nuclear target. The proton incident Drell-Yan reaction'^^l 
on nuclear targets provides, in particular, the possibility of probing the propagation of 
projectile through nuclear matter, with the produced lepton pair not interacting strongly 
with the partons in the nuclei. Therefore, the nuclear Drell-Yan process is an ideal tool to 
study initial state energy loss effect. At parton level, the nuclear Drell-Yan process is only 
sensitive to the quark energy loss'^^^. However, a charm quark pair is predominantly due to 
gluon fusion in proton induced J/%1) production on nuclear targets. The gluon energy loss 
can not yet be constrained by the experimental data though some theoretical papers'^^' on 
J /if) production employ the conclusion that the gluon energy loss is 9/4 times larger than 
the quark energy loss due to the difference in the color factors. As a result, the initial state 
energy loss effect is not pinned down in production which is predominantly due to gluon 
fusion. 

In the previous article P^l, by using the nuclear parton distribution functions from the 
global analysis, the nuclear Drell-Yan production cross section ratios were investigated for 
SOOGeV protons incident on a variety of nuclear targets in the framework of Glauber model. 
At hadron level, the energy loss of the beam proton in nuclear environment was determined 
well by fitting the nuclear Drell-Yan data from the Fermilab E866 experiment ^^^1. In this pa- 
per, by combining three representative sets of nuclear parton distribution functions with the 
initial state energy loss of the incident proton determined by the nuclear Drell-Yan reaction, 
J/ ^0 production cross section ratios for SOOGeV protons incident on tungsten and beryllium 
nuclear targets are analyzed with meanwhile taking account of the nuclear absorption effect 
in the Glauber model. It is hoped to have a new knowledge about the nuclear effects on 
J/ip production in proton- nucleus collisions. 

The paper is organized as follows. In sect.II, a brief formalism for the differential cross 
section in J/-^ production is presented. The section III is our results and discussion. The 
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summary is given in sect. IV. 



II. THE FORMALISM FOR J/ip PRODUCTION DIFFERENTIAL CROSS SEC- 
TION 

At leading order(LO) in pcrturbative QCD, the charmonium production cross section is 
the sum of two partonic contributions, gluon fusion [gg) and quark-anti-quark annihilation 
{qq)i convoluted with the parton distribution functions in the incident proton p and the 
target nucleus A^^l; 



da r"''^ , 2m 

dm 



f 



fg{xi,m ) fg(x^,m) agg(m ) 



q=u,d,s 

where Xi(^2) is the projectile proton ( target) parton momentum fractions, xp = Xi — X2, 
^/s is the center of mass energy of the hadronic collision, = X1X2S, = 1.2GeV and 
rrijy — 1.87GeV are respectively the charm-quark and D meson mass, and agg ( agg) is the 
LO cc partonic production cross section from the gluon fusion(quark-antiquark annihilation). 
Pj^^ is the fraction of cc pair which produces the J/ip state, ff and //^ stand respectively 
for the parton distribution function in the proton and in the nucleus. 

According to Glauber model'^^', the projectile proton scattering inelastically on nucleus 
(A) makes many collisions with nucleons bound in nuclei. The probability of having n 
coUisions can be expressed as 

n=l 



where 



'ml 



Cm(~ 30m6) is non-diffractive cross section for inelastic nucleon-nucleon collision, and T(o) 
is the thickness function of impact parameter h. 

Now let us take account of the energy loss of the projectile proton moving through the 
target nuclei before producing the charm quark pair in proton-nucleus. The energy loss due 
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to beam proton can induce the decrease of center of mass system energy of the nucleon- 
nucleon coUision producing cc, and affect the measured J/%1) production cross section. After 
the projectile proton has n coUisions with nucleons in nuclei, it is supposed for convenient 
calculation that the center-of-mass system energy of the nucleon-nucleon collision can be 
expressed as 

Vs- (n- 1)AV^, (3) 

where A-^/i is the center-of-mass system energy loss per collision in the initial state(see 
Ref.[21] for more detail discussion). Therefore, the J/ip production cross section in the nth 
collision can be rewritten as 

d(7(") _ da 
dx„ dx' 
Here the rescaled quantities are defined as 

x'p = rsXF, x[ 2 = rsXi,2, (5) 



^ (6) 



with the centre-of-mass system energy ratio: 



Furthermore, the nuclear absorption of the charm quark pair is added in final state effect. 
In the framework of the Glauber model, the probability for no interaction (or "survival 
probability") of the J/ijj meson with the target nucleus f^*^^ , can be calculated as 

= / dbil - e-(--)-(^)^r ). (7) 

The survival probability depends obviously on both the atomic mass number A of the target 
nucleus and the J/V^-nucleon inelastic cross section a'^lf . If the factorization between the 
cc production process and the subsequent possible J/ip inelastic interaction with nuclear 
matter is assumed, the J/ip production cross section in the nth collision is written as 

^ = —-Sabs. (8) 

dxp dx'^ * ^ 

Combining above ingredients on initial and final state effects, the J/ip production cross 
section in proton-nucleus collisions can be expressed as 

, da , , da^^"* , , 
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III. RESULTS AND DISCUSSION 



The Fermilab Experiment866(E866) W measured the differential cross section ratios of 
proton induced tungsten to beryUium target for J/%1) production, 

\ dx I 

Rw/Beixp) = , a„p-Be . - (10) 

The covered kinematical range was —0.1 < xp < 0.95. The experimental data were provided 
in small xp ( —0.1 < xp < 0.3), intermediate xp ( 0-2 < xp < 0.65), and large xp ( 0.3 < 
Xp < 0.95), respectively. By using three sets of leading order nuclear parton distribution 
functions together with CTEQ6L parton density in the protonl^^l, meanwhile taking account 
of the energy loss of the beam proton in initial state effect and the nuclear absorption of the 
charmonium states traversing the nuclear matter in the uniform framework of the Glauber 
model, a leading order phenomenological analysis is given in the color evaporation model 
to the E866 experimental data on the ratios of J/ijj production differential cross section 

Rw/Be{xp). 

As for the nuclear parton distributions in our calculation, the nuclear modification for 
gluon distribution function is apparently different between the different sets. The Figure 
1 shows the ratio of the tungsten nucleus over the proton parton distribution function at 
the charm quark pair mass scale, FC^{x,Q'^ — (2mc)^), as a function of Bjorken x for the 
leading order nDS (dotted line), EKS98 (solid line) and HKN07 (dashed line) nuclear gluon 
distribution functions. The main differences between the three nuclear paramctcrizations are 
found that nDS and EKS98 nuclear gluon densities have respectively nuclear anti-shadowing 
effect in the region 0.03 < x < 0.25 and 0.03 < x < 0.38. It is emphasized that nDS nuclear 
gluon density have comparatively very small anti-shadowing than EKS98 in intermediate x 
region, whereas HKN07 nuclear gluon distribution function has only the nuclear shadowing 
effect in small x. The ratio FC^ from HKN07 increases as x becomes larger, and crosses the 
line i?* = 1 at a; ~ 0.1. 

If neglecting the energy loss in initial state effect and the nuclear absorption in final 
state effect, the calculated results are compared with the E866 experimental data^^^ on J/ if) 
production cross section ratios Rw/Be{xp) in Fig.2. The solid, dotted and dashed fines 
are the results on Rw/Seixp) by using the EKS98, nDS and HKN07 parameterizations, 
respectively. As can be seen in Fig.2, RwjBeixp) from nDS nuclear parton distributions 
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FIG. 1: The ratio Rg of the gluon distribution in a tungsten nucleus over that in a free proton. 
The solid, dotted and dashed lines correspond to the results given by the EKS98, nDS and HKN07 
nuclear parton distributions, respectively. 

shows a remarkably flat behavior as a function xp in the region \xp\ < 0.1. However, 
Rw/Beixp) by HKN07 parameterizations turns into a tendency towards enhancement in the 
negative xp region. The anti-shadowing enhancement is appeared in the EKS98 model at 

~ 0. It is worth noting that the nuclear suppression from the nuclear effects on the 
parton distributions becomes larger as the increase oi xp in the range xp > 0. With taking 
example by HKN07 nuclear parton distribution functions, the suppression is approximately 
4% to 14% for Rw/Be{xp) in the ranges 0.1 < xp < 0.95. 

The curves in Fig. 3 show the calculated Rw/Be{xp) for J/ip production, compared with 
the E866 data, given by using the EKS98, nDS and HKN07 nuclear parton distributions 
together with the energy loss of the beam proton in initial state effects. In our calculations, 
the center-of-mass system energy loss per collision Ay^ = 0.18GeV was determined from 
the nuclear Drell-Yan experimental data in the Glauber model^^l . It is found that the energy 
loss effect of the beam proton in initial state has a small impact on the differential cross 
section ratios Rw/Be{xp) in the region xp < 0.2. However, the energy loss effect adds further 
the nuclear suppression on Rw/Be{xp) in the range xp > 0.2. The nuclear suppression on 
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FIG. 2: The J/'4> production cross section ratio Rw/Bei^p)- The fihed diamonds ( boxes and 
triangles) are the E866 experimental data W in the region —0.1 <xf < 0.3 ( 0.2 <xf < 0.65 and 
0.3 < < 0.95). With assuming only the nuclear effects on the parton distribution functions, 
the solid, dotted and dashed lines correspond to the results from the EKS98, nDS and HKN07 
parameter izations, respectively. 

Rw/Be^xp) from energy loss effect increases gradually in the region 0.2 < xp < 0.8, and 
becomes much steeper in the region xp > 0.8. The suppression due to energy loss effect is 
approximately 1% to 14% and 14% to 50% in the ranges 0.2 <xp< 0.8 and 0.8 <xp< 0.96, 
respectively. As for HKN07 nuclear parton distributions, the total suppression from the 
nuclear effects on parton distribution functions and energy loss effect is roughly 4% to 27% 
and 27% to 65% in the ranges 0.2 < xp < 0.8 and 0.8 < xp < 0.96, respectively. The similar 
results can be obtained from nDS and EKS98 nuclear parton distributions. Therefore, the 
energy loss effect, resulting in the suppression on Rw/Bei^p)-! is more important than the 
nuclear effects on parton distributions in high xp region. Although the initial state effects on 
Rw/Bei^p) become larger as the increase of xp, especially in high xp region, the remained 
deviation from the E866 data need to be contributed by the final state effect. 

After considering the initial state effects, it is supposed further that the nuclear absorp- 
tion of the charm quark pair is the main final state effect. The ratios of J/ if) production 
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FIG. 3: The ratio Rw/Bei^F) by combining the nuclear effects on the parton distributions with 
the energy loss effect in initial state. The comments are the same as Fig.2. 



TABLE I: The J/V'-nucleon inelastic cross section a^^f and x^/d-o.f. in the small xp region. 





HKN07 


nDS 


EKS98 




2.30 


2.01 


3.60 


xVd.o.f. 


4.371 


1.066 


11.758 



differential cross section Rw/Beixp) are calculated in the small Xp region by combining the 
nuclear effects on parton distributions with the energy loss and nuclear absorption effect in 
the uniform framework of the Glauber model. The theoretical results are shown in Fig.4 
against the E866 experimental data on Rw/Beixp)- The solid, dotted and dashed lines 
stand for the results from the EKS98, nDS and HKN07 parameterizations, respectively. 
The J/t/'-nucleon inelastic cross section a'^[f (in unit mb)and the P^r degrees of freedom 
are calculated and summarized in Table I by fitting the experimental data. 

As can be seen in Fig.4, the E866 experimental data points show remarkably flat behavior 
as a function of xp in small xp- The covered kinematical range —0.1 < xp < 0.3 correspond 
roughly to the parton momentum fraction 0.02 < x < 0.12 in target nuclei. The target 
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FIG. 4: The J/ip cross section ratio R{xf). The sohd circles are the ratios on tungsten to beryUium 
target from the E866[^l. The open circles are the ratios on tungsten to carbon target from the 
HERA-B experiment [^1. The solid, dotted and dashed lines are the calculated Rw/Bei^p) from the 
EKS98, nDS and HKN07 parameterizations by considering the initial state effects and final state 
nuclear absorption effect with a constant cr'^lf ■ 

parton momentum fraction increases as the decrease of xp in the negative xp region. In 
the range —0.1 < xp < 0.3, the gluon fusion channel dominates over the qq annihilation 
process in J/ijj production. By compared this observation with the predictions given by 
the different nuclear parton distribution functions, the tendency reported by the E866 data 
is well reproduced by the nDS and HKN07 parameterizations in small xp domain. These 
two parameterizations do not show a strong anti-shadowing effect, whereas the strong anti- 
shadowing presents in EKS98 nuclear gluon distribution at x ~ 0.1. It is shown that the 
E866 experimental data consequently tend to favor nDS and HKN07 sets rather than the 
EKS98 parameterizations. Recently, the HERA-B measurement indicates that the nuclear 
suppression on J/ip production cross section ratios turns into a slight tendency towards 
enhancement in the negative xp region. In Fig.4, the open circles show the ratios Rw/ci^F) 
(the Table 3 in Ref. [27]) which are derived from the HERA-B a{xp) measurement (the Table 
5.1 in Ref. [6]). It is apparent that the enhancement tend is only produced by HKN07 
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FIG. 5: The J/^-nuclcon inelastic cross section CJ"^^^ ci function of xp by means of HKN07(soli(i 
circles), EKS98(filled triangles) and nDS(empty circles) nuclear parton distribution functions. 



paramcterizations, which do not present any anti-shadowing effect. Therefore, it is concluded 
that the current experimental data is not in support of the anti-shadowing effect in nuclear 
gluon distribution. Our analysis shows similar tendency with the results in Ref.[26] in the 
anti-shadowing part by considering E866 measurement, but the new HERA-B data further 
rule out nDS nuclear gluon distribution with weak anti-shadowing effect. The papcr'^^l used 
the color evaporation model, and investigated J/ip production cross section ratio Rw/Be{x) 
as a function of the target momentum-fraction x(0.02 < x < 0.1) from ESGG'^^ by using 
the various nuclear parton distributions and keeping a'^^f as a free parameter. Their result 
without considering HERA-B measurement is in favor of the nuclear gluon distributions 
which do not predict a strong anti-shadowing at x ~ 0.1 nor a large shadowing at small x. 
In detail, the E866 data agree with nDS and HKM t^^l (without nuclear gluon anti-shadowing 
effect similar to HKN07)parametrizations, rather than the EKS98 nuclear gluon distribution. 
In addition, the a'^lf is slight bigger than our results in Table I because of not including the 
energy loss effect in initial state. 

The calculations on Rw/Bei^p) extended to the full experimental range inxp indicate that 
a single value of a'^lf does not provide a good agreement with the E866 experimental data. 
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Therefore, the J/^'-nucleon inelastic cross section is obtained as a function of by 
fitting the experimental data on a versus xp- In Fig.5, The sohd circles, filled triangles and 
empty circles correspond to the values of a'H'^ with based on three different nuclear parton 
distribution functions. In the region Q < xp < 0.2, the values of given by nDS are 
roughly a constant because of the flat nuclear modification on gluon distribution function. 
However, the calculated cr^^^^ from EKS98 have a steeper increase than those from HKN07 
as the decrease of xp, which originates from the nuclear modification on gluon distribution 
by EKS98 has sharply increase than that by HKN07 parameterizations in the corresponding 
region of the target parton momentum fraction. In the region xp > 0.2, the value of a'^[f 
comes to be larger as the increase of xp if excluding the two values with deviating from the 
tendency in the range xp > 0.7. With difference from the conclusion in Ref.[27], the a'^^f 
predicted here is smaller than that without the energy loss effect in initial state in high xp 
region. The results for a'^lf from HKN07 nuclear parton distributions can be represented 
for convenience by the simple parametrization: cr'^lfixp) — 2.685 — 9.038xf + 36.217x|,. 
Therefore, it can be confirmed that a'^lf depends on the kinematical variable xp. The 
neglect of energy loss effect in J/ip suppression studies is insufficient to gain an insight into 
the nuclear effects in J/ip production. 



IV. SUMMARY 

The precise identification of J/ip suppression mechanism as signatures of quark-gluon 
plasma formation requires a detailed and quantitative understanding of the nuclear effects 
already present in proton-nucleus collisions. The complex nuclear effects include mainly 
the nuclear effects of parton distribution functions and energy loss effect in initial state 
together with the nuclear absorption of the charm quark pair in final state. Although the 
nuclear effects on the parton distribution functions have been investigated by the global 
analysis analogous to those of the free proton, the nuclear gluon distributions have large 
differences between the different sets. By using the Glauber model, the energy loss of 
the beam proton in nuclear environment was determined well by fitting the nuclear Drell- 
Yan data from the Fermilab E866 experiment for SOOGeV protons incident on a variety 
of nuclear targets '^^l The nuclear absorption effect remained in final state can be studied 
with J/ip production in 800GeV protons incident on nuclear targets. We perform a leading 
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order phenomenological analysis on J/%1) production cross section ratios Rw/Be{^F)-i and 
compare with the E866 experimental data. The different nuclear effects are obtained on 
the J/ij) suppression. It is shown that the energy loss effect in J/tp suppression is more 
important than the nuclear effects on parton distributions in high xp region. In the small 
xp (—0.1 < Xp < 0.3) range, the E866 data are in agreement with nDS and HKN07 
rather than the EKS98 parameterizations which give a large anti-shadowing effect. However, 
the new HERA-B measurement do not support the anti-shadowing effect in nuclear gluon 
distribution, which further rules out nDS nuclear gluon distribution. It is found that the J/V'- 
nucleon inelastic cross section cr^^^ depends on the kinematical variable Xp^ and increases as 
Xp m the region xp > 0.2. Of course, the precise experimental data on J /ip production will 
be mandatory in proton- nucleus collisions, especially in the negative Xp region. Therefore, 
it is desirable to operate precise measurements at J-PARC^^^l and Fermilab EQOGf^^^ in the 
future. 
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